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■ Abstract 

■ Kinetic equations for ultrarelativistic particles with due account of gravitational inter- 
actions with massive particles in the Robertson- Walker universe are obtained. On the basis 

, of an exact solution of the kinetic equations thus obtained, a conclusion is made as to the 

■ high degree of the uniformity of the relict radiation on scales with are less than 10'. 



^ 1 Introduction 
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Bisnovatyi-Kogan and Shykhman (1982) have shown that within the scope of a macro-scopically 
homogeneous isotropic cosmological model based on the Newtonian theory of gravitation we 
can form kinetic equations with an integral of gravitational collisions which converge at large 
distance of action. The kernel of the collision integral obtained by Bisnovatyi-Kogan (1982) is 
the same as that of the Landau collision integral in Landau (1937). However the 'Coulomb' 
logarithm turns out to be finite at any finite cosmological time. Kandrup (1982) substantiated 
(f) . this result, but the collision integral differs from the one obtained in Bisnovatyi-Kogan and 

Shykhman (1982) in so far it does not vanish on account of equilibrium distribution. The 
difference between the kernels of collision integrals in Bisnovatyi-Kogan and Shykhman (1982) 
and Kandrup (1982) seems to be caused by the following. A Robertson- Walker universe lacks 
homogeneity in time, which is why the energy of the particles is not an integral of motion. On 
the other hand, in consequence of the long-range character of gravitational interaction the act 
of gravitational collision is protracted. It is precisely the combination of the temporal non- 
locality of gravitational interaction and the absence of temporal homogeneity that leads to the 
collision integral in Kandrup (1982). In Zakharov (1984) an integral of Coulomb collisions is 
obtained for nonrelativistic charged particles in a Robertson- Walker universe filled with dust. 
This integral coincides with the collision integral obtained by Kandrup, if the same notations are 
used, although all calculations in Zahkarov (1984) are made within the framework of Einstein's 
theory of gravitation. It should be noted that the employment of Einstein's theory of gravitation 
for the problem of a Coulomb interaction of non-relativistic particles is unwarranted, since 
non-relativistic particles interact by means of Coulomb's field in an identical manner both in 
Einstein's and in Newton's treatment, whereas the Coulomb component of an electromagnetic 
field in an isotropic space, being bound up with the law of the conser- vation of a charge, does 
not differ from the traditional component. The are precisely these facts that account for the 
agreement between the results in Zakharov (1984) and Kandrup (1982). 

Thus, the derivation of kinetic equations without regard to the radiation for non-relativistic 
particles in a Robertson- Walker universe can always be carried out within the scope of the 
Newtonian theory of gravitation and Newtonian mechanics. This assertion is based on two 
facts: (1) the adequacy of describing the effect of a gravitational field on non-relativistic particles 



1 Original article is published in [T|. In the given version original typing errors are corrected. 
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within the scope of Einstein's and Newton's theories of gravitation, which makes it possible to 
employ Newtonian mechanics; (2) the smallness of the path of a non-relativistic particle as 
compared to the light horizon, which makes it possible to employ the Newtonian theory of 
gravitation. 

There exists a fundamentally different condition in the case of ultrarelativistic particles. A 
descrition of the gravitational effect on these particles can only be carried out within the scope 
of a relativistic theory of gravitation. Indeed, on the one hand, the motion of an ultrarelativistic 
particle is affected by the components of a three-dimensional metric, which is why Newtonian 
mechanics cannot be employed to describe motion. On the other hand, an ultrarelativistic par- 
ticle runs a distance comparable to the light horizon, which makes the Newtonian description 
of a gravitational field unacceptable. However, in the case of electromagnetic interaction the 
question of forming kinetic equations for ultrarelativistic particles in a Robertson- Walker uni- 
verse is easily solved; in consequence of the conformal invariance of field equations and motion 
equations, these equations will be not different from those in a plane space (see Ignat'ev, 1982). 

This brings to the fore the basic problem of formulating kinetic equations for ultrarelativistic 
particles in a Robertson- Walker universe with allowance for their gravitational interaction with 
non-relativistic ones. Suitable techniques for forming these equations are provided by Ignat'ev 
(1978, 1983). In the present paper a kinetic equation for ultrarelativistic particles is derived, 
with due account of their gravitational interaction with massless ones. This is carried out on the 
basis of averaging a collisionless kinetic equation over the local fluctuations of the Robertson- 
Walker metric caused by the gravitational fields of massive point particles. The solution of 
Einstein's equations linearized about a Robertson- Walker solution shows that the local fields of 
point masses formed by a redistribution of Robertson- Walker matter, always remain within the 
sound horizon. In this connection the problem of the convergence of the collision integral at 
large distance of action becomes irrelevant. The kinetic equation thus obtained has a collision 
term whose structure is the same as that of the Belyaev-Budker collision integral (cf. Belayev 
and Budker, 1956), if one of the particles is taken to be non-relativistic in the latter. A solution 
of the kinetic equation thus obtained is presented. The solution describes the process of the 
isotropization of a homogeneous, but anisotropic distribution of massless relict particles. The 
estimates show that angular harmonics with a scale less than 10 angular minutes are strongly 
damped. All the notation, unless otherwise specified in the text, are according to those in 
Ignat'ev (1978, 1981, 1982, 1983), Belayev and Budker (1956), or Landau and Lifshitz (1972). 

2 Massive Particles in a Robertson- Walker Universe 

The gravitational field of a massive point particle in a Robertson- Walker universe is basically 
obtainable with the aid of the well-known Lifshitz solutions (see Landau and Lifshitz, 1972). 
However, the formulation of initial and boundary conditions in a problem with spherical sym- 
metry in terms of plane waves loses its physical lucidity. Besides, the synchronous frame of 
reference by Landau and Lifshitz (1972) is unsuitable for our purpose. We shall write the 
energy-momentum tensor of a massive point particle (Ignat'ev, 1983) as (h = = G = 1): 




(2.1) 



where D(x/x') is the invariant four-dimensional function of Dirac: i.e. 




2 



(dfl = yj—gd 4 x); integration in (2.1) is carried out along the entire world line of the particle 
x' = x'(S). As will be seen from the following, in a non-empty space the mass of a particle 
cannot remain constant; therefore, m = m(x) is a scalar function. The motion equations of a 
point particle with a variable mass is obtainable from an invariant Hamiltonian (Ignat'ev, 1982) 

H m (x,P) = —g ik P i P> e -m = 0, 
m 

from which we find using standard procedure (as in Ignat'ev, 1983) 
du l 



dS 



= (lnm), fc (g ,fc -wV). 



(2.2) 



If we assume the spherical symmetry of the problem we shall write the space-time metric in 
isotropic coordinates 



dS 2 



e v dv 2 - e x [dr 2 + r 2 {d6 2 + sin 2 9dip% 



(2.3) 



where v = v{r 1 rj), A = A = (r, r/). Then the world line of the particle is the time line r — 0, m 
m{r}) being an arbitrary function in consequence of (2.2). Integrating (2.1) we get 



5Tl = 6\6ime^ 2X 8{v) 



where S(r) is Dirac's three-dimensional S - function in plane space. We shall represent the total 
energy-momentum tensor T lk as 

rjiik rpik _|_ fij^ik 

where Tj k is the energy-momentum tensor of a fluid. In consequence of the Einstein equations, 
(2.1) and (2.2) the law of conservation is satisfied by 

Tf k = -V 4 J mD(x\x')dS. 

In particular, it is, in our case, of the form 

Tf k = -me-(" +A '%). 

Einstein's non-trivial equations for the metric (2.3) take the form (see Landau and Lifshitz, 
1972) 



1 

r 



^ + AV + -(A' + v') 



1 • 



3, 



e- v [X--\v+-\ 2 ) =87r\p + v 2 (e+p)]; 



—e 



2{\" + v") + v !2 + -{\' + V ') 
r 



8tt 



me~i x S(r) + e + v 2 (e + P) 



V (2A' - u'X) = 8ir(e + p)e^ v+X) vV~l + v 2 



where v = u r e x ' 2 is the reference projection of the radial velocity of a fluid. If we subtract the 
second equation from the first, we get 



r 



^- + AV-^- + i(A' + z/)-(A" + */') 



87r(e + p)v 2 . 
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We shall take m, v, \' ,v' to be first-order infinitesimals. Then in a linear approximation the 
last equation is easily integrable, yielding 

\ + v = C 1 (r 1 )r 2 +C 2 ( V ). 

By making a direct substitution we can very that C\ = 0. Admissible transformations of 
coordinates which preserve the form of the metric (2.3) are 

rj = 7](fj); r = kf; (k = Const). 

In consequence, we can add to v an arbitrary time function, and to A an arbitrary constant. 
Let us select this function so that 

A = In a 2 +£(r, n), v = In a 2 - £(r, 77), (2.4) 

where a = a{n), £ <C 1. Linearizing Einstein's equations with respect to the smallness of m, v, £, 
we obtain the system 

1 (2a d 2 \ > 2 „ 

_(_-_j=- 87 rp ; 3^=8^; (2.5) 

i + 3i* + (™-*)t = -***pe-, (2.6) 
a \ a a z J de 

' 9 (r 2 a + 341-K) - 8^ + 8^fe; (2.7) 



a 2 r 2 Q r a 4 Q n a 3 

1 d 



8ir(e +Po) drj 



«), (2-8) 



where po = Po(f?) ; £ o = £o(j])> an d we have put 5p = (dp/dede. Equations (2.5) describe the 
evolution of a flat-space Robertson- Walker universe (cf. Landau and Lifshitz, 1972); Equation 
(2.8) determines the radial velocity of a fluid. To solve the singular equations (2.6) and (2.7), 
let us put 

£=— (m-*), (2.9) 
ra 

and m = m(r)), 

lim <J(r, r))/r < +00. (2.10) 

r— >0 

If we substitute (2.9) into (2.7) and taking account of (2.10) we get 

4ira 3 6e = 1*" + — A( m _*). (2.11) 
By use of (2.11) in the right-hand side of Equation (2.6), we get closed equations for m{rj) and 

a ■ ( ,Av\ /2a d 2 \ 

* + £*fi + 3^+*f--4 > )=*"?- ( 2 - 13 ) 

a \ de I \ a a A I de 
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At the non-relativistic stage p = and, according to (2.5), a ~ if. Then (2.12) and (2.13) are 
easily integrable: yielding 

m = a(^)\,{£); * = W {r)(*)\v{r){£) , (2.14) 

where a, \i are arbitrary constants and W(r),V(r) are arbitrary functions. Thus, the mass of 
a point particle in a medium is not conserved, which is the consequence of its gravitational 
interaction with the fluid. 

To determine the functions W(r),V(r) it is necessary to consider the evolution of a metric 
at the ultrarelativistic stage of expansion (e = 3po)- In this case Equations (2.12) and (2.13) 
take the form 

2 2 2 • 2 1 

m+-m 7rm = 0; * + = . 

rj r] z r\ r\ A 3 

The first one is easily integrated, and yields 

ro=0 /f*y+^ i (2.i5) 

V V J Va 

where c',/z' are arbitrary constants; to solve the second equation we shall make a substitution 
r\ = V3t and 

*=^--$(r,T). (2.16) 
or r 

On substituting (2.16) into the second equation and changing the order of differentiation we 
shall reduce it to the form 

Ai($ rr -$ rr ), (2.17) 

OT T 

whence 

$(r,r) = $ + (r + r) + $_(r,r). (2.18) 

To the solution obtained we can add rF(r), where F(r) is an arbitrary function. However, this 
addition, according to (2.16), does not alter the value of *ff. A particular solution to (2.18) 
is $ = 0. Then, according to (2.9), we shall obtain a Newtonian potential caused by a point 
particle of variable mass m(rf) determined with the aid of (2.15). However, this solution is 
physically unacceptable for the description of the gravitational field of a point particle which 
arose as a result of fluctuation at the moment of time t — 0, since it is inconsistent with 
the principle of causality. Indeed, the Newtonian potential referred to make it possible to 
obtain information about a particle beyond its light horizon and even to determine the mass 
of the particle. Therefore, boundary conditions should be formulated in such a way that at 
least beyond the light horizon the potential £ should vanish together with its derivatives. Such 
boundary conditions are in according with the 'birth' of a particle as a result of the redistribution 
of Robertson- Walker matter. In fact, however, the horizon of the potential is not the light cone 
but the sound cone, since the change in the potential is bound up with the redistribution of 
matter, and the latter proceeds precisely at the speed of sound. Mathematically this is related 
to the characteristics of Equation (2.17) and it is precisely at the sound horizon r = t = r/y/3 
that boundary conditions take the simplest form. 
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For the mass of a particle to remain limited at t — > it is necessary to put a' — in the 
solution (2.15). Boundary conditions at the light horizon which satisfy the principle of causality 
are of the form 

£(r,T)| r=r = 0; C(r,T)\ r=T =Q; v(r, r)| r=r = 0. 

The latter condition is, however, automatically satisfied in consequence of the first two condi- 
tions. If we substitute £ here from (2.9) we shall reduce the first two conditions to the form 



V{T,T)=m(T)=LL-- *(r,r)| r=T =0. (2.19) 

TO 

According to the definitions of (2.17) and (2.18) these conditions are equivalent to 

$+(2r) + $_(0) = 2^^T 3 + Ar; (2.20) 

TO 

&+(r + r)| r=T - $'_(r - r)| r=r = Br, (2.21) 

where A and B are arbitrary constants. Besides, due to (2.10) another condition is to be 
satisfied: namely 

$+(t) + $_(t) =0. (2.22) 
From (2.20) and (2.22) we find that 

$ + (x) = ~^x 3 + ^x + <£ + (0), 

4 T I 

$-(x) = -~^x 3 --x-<S> + (0). 
4 t 2 

Differentiating these relations and substituting the results into (2.21), we get A = 0, = — ////Yq. 
If we substitute the obtained value of (, ) into (2.16), we finally gecl 




r < r; 

¥(r,T) = { * ''° \ (2.23) 

r > t. 

Let r]o be a moment of 'time' when the ultrarclativistic stage is replaced by the non-relativistic. 
Suppressing in the solution of (2.14) terms which correspond to the dispersion of the mass and 
joining this solution to (2.23) at the moment rj 0l we get an expression for the potential £(r, r) 
at the non-relativistic stage of expansion 



2fi 

e(r,r)=< ~ 



r ( r 2 
2t V t z 



r < t: 

(2.24) 



0, r > t. 

The potential obtained is time-independent. Therefore, if the distance Ar from a certain 
observer, synchronous in a Robertson- Walker metric, to a massive particle is more than the 
value r]o/V3, the observer will never be affected by the gravitation of a particle. Thus, in a 



2 A solution to (2.23) is also obtainable as an auto-model one, putting \t = r\J/(r/r). 
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Robertson- Walker universe the effective range of the gravitational forces of point particles turns 
out to be finite. The ratio between this range l g — a(rj)rj /V^ and the distance to the light 
horizon l c at the non-relativistic stage of expansion is less than unity and decreases as time goes 
on 



lc VS\t 



(2.25) 



At to ~ 10 13 c, t ~ 10 18 c, lg/l c ~ 10~ 2 . It is of interest to note that the dimension of gravita- 
tionally bound regions l g in the present epoch turns out to be of the order of 100 Mpc, which 
at the medium density g = 10~ 29 g • cm' 3 Mpc -1 indicates a mass of the order of 10 18 Af Q . It 
is also worth of noting that the mass m ~ m p i per a Planck moment of time increases to values 
of the order of 10 26 g in the present epoch. 

If we substitute Equation (2.24) into (2.11), we shall obtain an expression for the perturba- 
tion of the energy density caused by a massive particle, of the form 

- = -*^+«r)-^+«r), (2.26) 

from which follows an obvious fact that the perturbation of the energy density remains small 
as long as the mass of pa particle is small as compared to the mass of the gravitationally bound 
region: 

m(v) < 3^o« 3£ o(»?) = -| = m. (2.27) 

At to = 10 13 s we have m = 10 18 M Q . According to (2.26), at the boundary of the gravitation- 
ally bound region (r = rO) there arises an abrupt change in density 6e(r0, T))/eo(r)) — —m^/rfi, 
which vanishes at r\ — > and increases as time goes on. At the end the entire matter from within 
the sphere with the radius ro accretes by a massive partickll and the finite gravitational field 
in an empty sphere will be described by a Schwarzschild metric with the gravitational mass m. 
It is not difficult to obtain this finite metric by joining a Schwarzschild metric in an isotropic 
coordinate system to a Robertson- Walker metric (see, for instance, Landau and Lifshitz, 1972) 
on a sphere with the 'radius' ro = const. A smooth joining is possible precisely at r g = 2m, 
where m is described by Equation (2.27). 



3 Local Fluctuations of the Metric and Their Averages 

Let there be in a Robertson- Walker universe not one, but several massive identical particle^] 
with the coordinates v a — {x a , y a , z a }. Then the total space-time metric can be approximately 
written in a form (linear approximation with respect to m): 

dS 2 = {g lk + h lk )dx l dx k , (3.1) 

where go, is a Robertson- Walker metric, 

^fe = -a 2 ^^ea(|r-r a |). (3.2) 



3 If this is not impeded by the forces of pressure or a large angular momentum of the dust. 
4 If the masses of the particles are different, it is also necessary to average out the formulae obtained, with 
respect to the distribution of the masses. 
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Indeed, a contribution to the metric caused by the interaction of particles is of the order m 2 ; 
non-diagonal components of the metric tensor g a 4 ~ £v Q ~ m 2 are also of the same order. 
Introduce a field of observers, macroscopic in the metric of (3.1), whose coordinate grid is on 
massive particles. Such observers are geodesic with respect to the Robertson- Walker metric. 
Let the coordinates of massive particles r a assume random equiprobable values throughout the 
entire three-dimensional space, with no correlation between the positions of these particles. Let 
TV = Const, be the number of massive particles per Volume' V = 4/37rrg. Let us introduce an 
operation of averaging a certain field value ip(x\ri,r2, ■ ■ ■), which is a function of the positions 
of the massive particles, on the scale of macroscopic observers: i.e., 

(</?(#)) = IItF / d 3 T a (p(x\xi,x 2 ,...) (3.3) 

a Va J 

where integration is carried out within the spheres with the radius ro with centres at the point; 
V a — T^H Then the average of the local fluctuations of the Robertson- Walker metric of (3.2), 
according to (2.24) and (3.3), is 

(h lk ) = -a 2 S lk ^- = a 2 (r/)Const. (3.4) 
5r 

If we follow the procedure of Ignat'ev (1978), we shall renormalize the macroscopic Robertson- 
Walker metric and the local fluctuations hik so that the average of the latter should be equal 
to zero: i.e., 

Qik -> 9ik + (htk); hik -> h lk - (hik); £ Q -> - (3-5) 

5r 

Due to (3.4) a renormaliztion of the macroscopic metric gik reduces to multiplying (744 and g a p 
by constant numbers. Carrying out a further admissible infinitesimal scale transformation of 
the coordinates r\ and r with a constant scale coefficient, we shall restore the former value of the 
Robertson- Walker metric 

(dS 2 ) = a 2 (f])(drf - dx 2 - dy 2 - dz 2 ). (3.6) 

By use of the value of the function of £(r) redefined according to (3.5) we shall evaluate the 
change in the total mass SM(rO, 77) within the radius ro caused by a massive particle. For 
this purpose, we shall transform the expression for the total mass M(r) in the Schwarzschild 
coordinate system (Landau and Lifshitz, 1972), linearize it with respect to £ and employ the 
formula (2.24). As a result we get 

ro 

5M(r , rj) = Air— — / £,r 2 dr = 
a drj J 


Thus, the solution we have obtained actually describes the gravitational field of a particle formed 
by a redistribution of Robertson- Walker matter. 

If we calculate (£ 2 ( r )) according to (3.3) and (3.5), we get 



'•n 



(r(r)>=^(^(r)> = — N[^) . (3.7) 



5 In the formula (3.3) integration could be carried out over the entire space; however, this would yield nothing 
new, since according to (2.24) outside the sphere with the radius correlations vanish. 
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For the local and medium perturbations of the metric to be small, two conditions are to be 
satisfied: namely, 

\r a - r\ » 2/i; (3.8) 

7V (^) 2<<1 ' (3 ' 9) 

Since several massive particles can now find themselves within the sphere with the radius ro, 
the condition of (2.27) must be replaced by 

ft \ 2/3 

nm(r]) = Nfj, f — J < to. (3.10) 
By use of the properties of isotropy and homogeneity we can show the validity of the equalities 

{hikhlm) = SikSlmd 4 ^ 2 ); 

(3.11) 

(dj(h ik )hi m } = 0; (d4{hik)dj(hi m )) = 0. 

Averages of the type {d a hikdphi m ) diverge as r^ 1 , the divergence of these values is connected 
with the ordinary divergence of the particle's energy. To calculate these values we shall note 
that, due to the isotropy of space, 

Thus, 

((0aO(fyO> = 4^>- (3A2) 

By use of Equation (2.7) and the explicit form of the function A£ in (2.9), (2.24) to calculate 
£(r), we find 

-A£_ a = 8n f iS(r - r„) - ^. 



As a result we get 

-<^ 77 = ^ Js(r)Ur)dV. 



This integral diverges as r^ 1 , which is caused by the divergence of the total energy of the 
particle. To calculate the integral we shall employ the procedure of renormalizing the mass, 
treating S(r)/r as —S'(r). Then, taking into account (3.12), we get eventually 

6ttN (2fi\ 2 



<(aa0(^0> = -5-hr s °f>- 
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If we calculate the averages of corrections to the Einstein tensor caused by the fluctuations of 
the metric we shall find the corrections to the energy-momentum tensor of Robertson- Walker 
dust caused by the energy of the local gravitational fields: 

Thus, we shall obtain corrections to the energy density Se g and pressure Sp g of the forn[^ 

& ' = -^ = -4^(^) a ^ & ' + * = °- (3 ' 13) 

The allowance for the local fluctuations of the gravitational field in the Einstein equations is 
equivalent to adding to the Robertson- Walker dust a fluid with the equation of state e = —p. 
Note that the ratio Se g /so increases in proportion to rj 2 , but remains small when the conditions 
of (3.8)- (3. 10) are satisfied. 



4 Derivation of the Kinetic Equation 

Let f(x, P\x\, X2, ■ • •) = f(x, P) be a macroscopic function of the distribution of massless parti- 
cles; then the total number of particles recorded by the observers, associated with the velocity 
field u % on hypersurface S orthogonal to this field (cf. Ignat'ev, 1983) is 



/ d 4 Pu l P t 5(H)f(x,P), (4.1) 

where d A P = dPidP^dP-^dP^j ^J— g; H(x,P) is the invariant Hamiltonian of massless particles. 
We have used a bar to note the fact that all the values determined in a Robertson- Walker 
universe perturbed by local fluctuations. In accordance with our procedure we shall select 
as u 1 a field of observers which is geodesic in a macroscopic Robertson- Walker universe, then 
u 1 — SX/y/gZi ■ Thus, taking into account (3.1), (3.2), (3.4), and (3.5), we get from (4.1) 

L = III 1^1 dCipP * 5 ^^ P )' ( 4 ' 2 ) 
We shall take into account the fact that the Hamiltonian of massless particles is of the form 

= (4-3) 
where P 2 = P 2 + P 2 + P 2 . Transform the formulae (4.2) and (4.3) to the new variables x % , P,: 



p a = P a ; p * = Wrr|; ( 4 - 4 ) 



then 



H(x,P) = T ^H(x,P), (4.5) 



"In original article 1 because of negligence the error in a sign on density of energy Se g is admitted. However, 
expression for 8T^ is resulted the correct. 
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where H(x,P) is the Hamiltonian of massless particles in an unperturbed Robertson- Walker 
universe (3.6). Modifying the formula (4.2) taking into account the properties of 5(H), we get 

dh 



-mi 



d A P P 4 5(H)f(x,P), 



(4.6) 



which is precisely the same expression as the one in terms of an unperturbed metric (3.6). It 
should be noted that in (4.6) the fluctuation of the metric is present only in the function of the 
distribution of f(,P). Averaging (4.6) according to (3.3) we find that 



(L) 



IIH 



d A P u l P l 5(H)f(x,P), 



(4.7) 



where 



f(x,P) = (f(x,P)). 



(4.8) 



Thus, the average number of massless particles is evaluated with the aid of an ordinary formula, 
in a macroscopic Robertson- Walker universe, with respect to the distribution function averaged 
over the local fluctuations of the metric. This is the advantage of the frame of reference that we 
have selected. If we had selected, for instance, a macroscopically synchronous frame of reference, 
then instead of the simple relationship (4.8) we would have had a complex relationship which 
would include the derivatives of the distribution function f(x, P) and the correlations of this 
function with the local fields. This is caused by the fact that a macroscopic ally local observer 
is subjected to the effect of the local gravitational fields, as a result of which the scales of clocks 
and rods on the micro - and macro - levels are different. 

Now we shall proceed to average the microscopic collisionless kinetic equation (as in Ignafev, 
1981) 



5(H) 



dH Of dH df 



dPi dx l dx l dR 



= 



(4.9) 



If we carry out in (4.9) a preliminary transformation to the new variables x l ,Pi according to 
(4.4), we get 



5(H) 




dH_ df_ 

dP a 4 <9P 4 



dH_ df_ 

dP 4 A dP a 



(4.10) 



In Equation (4.10) all the momentum variables P a and P4 are treated as independent. Express- 
ing P4 in terms of P a with the aid of the mass shell equation H(x, P) — 0, we shall reduce 
(4.10) to a simpler form 




z a/ 1 / i + e 
t p *aj-2{ da ln w 



p 4 p i 



df_ 

dx a 







(4.11) 



(P a — —P a /a 2 ;P 4 = P4/ a 2 ). Our task is to obtain a kinetic equation for the macroscopic 
distribution function f(x, P) to within an accuracy of terms quadratic with respect to the local 
fluctuations of the metric. Therefore, we shall expand (4.11) into a series with respect to the 
smallness of £ confining ourselves to the second-order terms: 



dx° 



drj 



dx° 



0. 



(4.12) 
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Taking into account the fact that correlations between the positions of massive particles are 
absent, we shall represent the macroscopic distribution function f(x, P) in the form 

f(x,P) = f(x,P)+g(x,P), 

where 

g(x,P\x l ,x2, ...) ) (4-13) 

a 

and, in accordance with the definition (4.8), 

(g(x,P))=0. (4.14) 

If we substitute (4.13) into (4.12) and averaging the equation obtained, taking into account 
(4.14), we get 

Averaging (4.12) over all the particles except the ath one, taking into account the results of 
Equation (4.15) in the equation obtained, and suppressing terms quadratic with respect to £, 
we shall find an equation to determine the correlation function 

g a (x,P\x a ) = g a , (4.16) 

The integrals of Equation (4.16) are the integrals of geodesic lines in a Robertson- Walker uni- 
verse: namely, 

P a = Const; / = x% + n a {r)-r) ), (4.18) 

where n a = p a /p 4 = —P a /P 4 = Const. If wc integrate (4.16) along the trajectories of (4.17), 
we find that 

v 

9a = ~J dv'U\ r'(v')-r a \)-^jf[v'(r 1 ),r 1 ',P i \ + 

na 
n 

+P 4 J dr ] 'd a U\r'(v')-r a \)-^f[r'( V ),v',Pi]+ga(x a -7r a (r 1 -r 10 ),P i ). (4.19) 

We shall expand in the integrands of (4.18) the function f[r' (r]'), r]', Pi] 
/[rVW,^] - f(r,V,Pi) (V-V')y 4 §t. 

But according to (4.15) P l (df /dx l ) = 0(£ 2 ); therefore, the distribution function to within an 
accuracy here required, can be factored outside the integral in (4.18). For the same reason 
the function g a can be represented in the form of an arbitrary linear operator acting on f(x,P 
selecting it so as to satisfy the condition of (4.14). Having noted this, we shall proceed to 
calculate the integrals in (4.18). A non-zero contribution to these integrals is only provided by 
the regions 

I rV) -r„| = I r-r + 7?fa-»/)| <r„. (4.20) 
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A non-zero contribution to the averages in Equation (4.15) is provided by the regions 

|r-r a |s:ro. (4.21) 

Therefore, we are interested in the values of correlation functions at the intersection of the 
regions (4.19) and (4.20). 

To simplify the calculations we can put r a — and direct the velocity of the particle along 
the z-axis. The value of the correlation function obtained by integration is, generally speaking, 
time-depended; however, in the region (4.21) 

7] - Tfy > 2r , (4.22) 



Q 



2fi 



no 



9 
5r 



Q 



9r 3 
z 'o 



r a + z a 











^2)3/2 







* a = J d^drj' = £TT a + 2^(S af) - 7r Q 7r' 3 ) x 



no 



g 2 r a 



g 2 r 



It is necessary to generalize the expressions obtained, in the manner 
z a -> (7?,r - r Q ); g 2 = |r - r a \ 2 - (7? , r - r a ) 2 ; 
z a y |r r a |; x^ y x^ x^. 

Since (\I> a ) = 0, the correlation function can be represented as 

ga = -V {Q - m+ P<°L*~ 

Calculating the averages in accordance with (4.23), we get 

2 



(2 ± 



K: 



3N (2(i 



2r \ r o 



(6 a p - 7r Q 7r / g)A, 



where 



K 



17 



160 144 



7 n „ 29 
In 2 



31 



13 ^-^1-0.11. 



50 240 180 V 25 21 



(4.23) 



(4.24) 



(4.25) 



The last value of (4.24) is logarithmically divergent at the lower limit r —> 0, which is caused 
by the coincident illegitimacy of a linear approximation of the Einstein equations and a Born 
approximation. According to (3.8) at r ~ 2/j, the local fluctuations of the metric become 
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large, and at the same time the deflection angle of massless particle becomes large; therefore, 
integration in (4.24) is 



A 



2/x 



In2 + -. 



Substitute (4.23), with the acount of (4.24) into Equation (4.15) 



' drj dP a 
2r \r 



dP a 



APaP 4 - 



dP a 



Pa drj 
Pi I S a R - 



P a Pf, 

r 4 



df 



dPp 



= 0. 



If we change the order of differentiation in the fourth and fifth term, respectively, suppressing 
the small term which is proportional to ro/rj and transposing the final term to the right-hand 
side, we eventually obtain kinetic equation for massless particles in the region (4.21) 



a/ 

drj 



<df_ 
dx a 



3 N 
2^ 



r J dP a 



Pi 6. 



'a/3 



P a Pp 
p2 



df_ 

dPp 



.(4.26) 



The term on the right-hand side of Equation (4.25) is a collision 'integral', it describes the 
process of altering the momentum of massless particles by interaction with the local gravitational 
fields. The supplementary term on the left-hand side of (4.25), §(£ 2 ) P^d f / 'drj describes the 
effective change in the velocity of a particle in local gravitational fields. Indeed, if we suppress 
the collision term in (4.25), then the equations of characteristics will be of the form 



P a = Const; 



x = 



'-o ■ 



These equations describe the motion of a particle with the velocity 
5 



c 1 



< c. 



(4.27) 



A local observer who measures the velocities of massless particles equal to the speed of light and 
assumes that medium velocity of these particles coincides with the speed of light will observe 
a local violation of the law of conservation of the particle number. Indeed, if the observer 
knows that a massless particle escapes from a certain point in space in his direction he has 
every reason to expect to register it within a completely definite time At — Al/ c. If this fails to 
occur, the observer will record a local violation of the law of conservation of the particle number. 
However, since all massless particles travel at identical velocities (4.26), the observer can define 
this velocity as the speed of light. In agreement with this operation is the renormahzation of 
the momentum in (4.25), in which the term on the left-hand side of the equation, quadratic 
with respect to £, vanishes. 



5 The Laws of Conservation and Generalized Kinetic Equa- 
tions 

Carrying out the foregoing renormahzation of the momentum and changing over from a syn- 
chronous frame of reference, to an arbitrary one, we shall write Equation (4.25) in an invariant 
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form 



3 » fVi\„, a/ 



where 

2 r a V »"o / oPk 

W lk = (u,P)[PiP k +g lk {u 1 Pf - {u,P){uiPk-u k P^\. (5.3) 

The coefficient in (5.2) can be given a clear physical meaning by taking account of the determi- 
nation of the particle mass (to = /ia) and the particle number TV. Let us introduce the medium 
mass density g m = A r TO/4/37r(ar ) 3 , which corresponds to massive particles. It should be noted 
that this value changes as g m ~ i~ 4 / 3 at the linear stage of accretion as distinct from the total 
mass density g <~ t~ 2 . Then (5.2) takes the form 

I l = SnKg m W lk ^L. (5.4) 

The symmetric tensor Wi k has the following properties: 

W lk P k = 0; W lk u k = 0. (5.5) 
Let us integrate (5.1) over the space of momenta 



where 



J P\f 5{H)d A P 



is the density vector of the particle-number flux. Thus the local law of the conservation of 
particle number is also satisfied for a non-renormalized equation (4.25), which can be verified 
by integrating (4.25) over the space of momenta and the entire three-dimensional space. We 
shall multiply (5.1) by P l and integrate over the space of momenta. If we integrate twice by 
parts on the right-hand side of the equation, we get, taking into account (5.4) 

If = 32Ag m m(g* k - u'u^TW, (5.6) 

where 



T ik = J pipk j £(#) rf 4p 



is the energy-momentum tensor of massless particles. The energy of massless particles is con- 
served, which can be verified by comparing (5.5) with the time-like vector of conformal motion 

xi l = S\: 

^ = 0, 

where V 1 = T lk £ k . However, the three-dimensional momentum of the particle is not conserved, 
which is the consequence of the non-conservative nature of the system of massless particles, - a 
part of the momentum of massless particles is transferred to massive ones. 
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The kinetic equation (5.1) admits a natural generalization in the case where massive particles 
travel at arbitrary, but small velocities at which the rate of accretion is not altered significantly. 
Let F(r], P') be a function of the distribution of massive particles with masses which are now 
arbitrary, so that 



P')P' i d i P' 



is the density vector of the flux of massive particles. Then it is known that collision term of 
(5.2) (P- = raui) can be written as 

U = Snj Am 2 F(P')W lk ^d i P'. (5.7) 

To make sure of this it is sufficient to put F(P') = m~ 1 5(P)8(H m /2), i.e., to consider the 
distribution of fixed-mass particles at rest. The collision term of (5.6) takes account of the 
processes of transferring a momentum from massless particles to massive ones, but fails to take 
account of the reverse processes. To take into account the latter it is necessary to add to (5.6) 
a term which is antisymmetric respect to the transposition of particles 



L = 8tt / d i P'Am 2 W, k 



F(p') d 4^-f(p') dnn 



(5.8) 



dP k Jy ' 8P' k 

The kernel Wi k of the collision integral thus obtained is the same as that of the Belayev-Budkcr 
(1956) collision integral, if in the latter we let the momentum of one of the particles tend to 
infinity. The difference lies in the mutiplier 87r (instead of 2tt in Belayev and Budkcr (1956)), 
which is the consequence of the well-known effect, viz., the deflection angle of a photon in a 
gravitational field with the scope of Einstein's theory is two times as great as its Newtonian 
value. 

The collision integral of (5.7) as distinct from (5.1) now satisfies all the necessary laws of 
conservation. It can be used as the basis for showing that the total of energy and momentum 
within the system of 'dust + massive particles + massless particles' is conserved. 



6 Isotropization of Homogeneous Distribution of Massless 
Particles 

An exact solution of the kinetic equation (5.1) (as well as (4.25)) is an arbitrary isotropic 
distribution f(P). If the initial distribution is anisotropic, however, it will be isotropized by 
gravitational interactions. Let the distribution of massless particles at the momentum of time 
i] = rji in a synchronous frame of reference be of the form 



f( m ,P a ) = J2flm(Vl,P)Y l m (0,<fi), 



(6.1) 



l.m 



where 9 and <p are the azimuthal and the polar angles in the momentum space, respectively. 
Representing the distribution of f(r],P a )) in a form analogous to (6.1), and on separating the 
variables, we obtain equations for the functions fi m (r], P), 

dflm 



drj 



-1(1 + l)mAafi m g m 8Tr, 



flmiv, P) = fim(m,P) exp 



-1(1 + 1)8tt / Amg m dt' 



(6.2) 



16 



In particular, for the dependence of m(t) and g m (t) we have obtained, we shall find by inte- 
grating in (6.2): 



flm(V, p ) = flm(Vi, p ) ex P 



-1(1 + l)247rAm (?m [ 1 - ( — 



(6.3) 



where m = m(t) , g m = Q m (t) . At t — > oo the expression in the exponent increases proportional 
to t 1 / 3 . Therefore, at t — > oo the expression (6.1) retains only one harmonic with I — 0, i.e., the 
distribution is isotropized. The ^-harmonic is presented by the angular scale Aip = 2n(l + 1). 
Therefore, at the given present-time values the particle masses m(t) and their medium density 
Q m (t) all harmonics with the angular dimension 



Aip2n < y / 24:irAmg m t 



(6.4) 



vanish. Atm~ 10 16 M Q , g - 10 30 g cm' 3 and t = 2 x 10 10 yr we obtain from (6.4) Aip < 10 
angular minutes (Figures 1 and 2). 




Figure: 1. The dependence of the damping coefficient of the an- 
gular harmonic distribution 7 = fi m (to)/fim(t) from the angular 
scale A(p = 2tt(1 + 1). Curve 1: m = 10 16 M Q , g m /g = 0.1. Curve 
2: m = 10 16 M Q , g m /g = 0.2. The angle Aip is measured in 
angular minutes. 

Thus, any relict radiation must be highly homogeneous on scales less than 10 angular min- 
utes. This effect can resolve the contradiction between the deductions of a standard adiabatic 
theory of galaxy formation (as in Zeldovich, 1983) and observation data which give evidence in 
favour of the absence of small-scale fluctuations of relict radiation. 

Now let us assume that the accretion is completed at the moment of time ti i.e., m = Const, 
and g m <~ t~ 2 . Then in place of (6.3) we obtain from (6.2), 



flm(V, P ) = flm(m,P) exp 



t 2 

-1(1 + l)8TTAmg m — 
ti 



1 - — 



(6.5) 
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In case, at t ~ oo the exponential index becomes large, 4Al(l + l)m/3ti. This implies that har- 
monics with an angular scale of the order of tens of degrees can be dampened, and isotropization 
takes place according to (6.5) at the earliest stages following the termination of accretion. Ap- 
parently, this cannot be the case. But it follows then that the masses of super aggregations 
should be either less than 1O 16 M , or their formation should be completed at significantly later 
moments of time t\ > 10 16 s. 

A^P (angular minutes) 



5 





2 




O 


? ■ 


'55 


5- 


dimei 


2 ■ 


to 
u 
c 




m 



05 1 1,5 Z 

Figure: 2. The dependence of the maximum resolution of Atp on 
the red shift Z: m = 10 15 M Q ; g m /g = 1. The graph shows the 
values of Af for which 7 = . 
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